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Zurek suggested (Nature 317, 505; 1985) that the Kibble mechanism, through 
which topological defects such as cosmic strings are believed to have been cre- 
ated in the early Universe, can also result in the formation of topological de- 
fects in liquid 4 He, i.e. quantised vortices, during rapid quenches through the 
superfluid transition. Preliminary experiments (Hendry et al, Nature 368, 
315; 1994) seemed to support this idea in that the quenches produced the pre- 
dicted high vortex-densities. The present paper describes a new experiment 
incorporating a redesigned expansion cell that minimises vortex creation aris- 
ing from conventional hydrodynamic flow. The post-quench line- densities of 
vorticity produced by the new cell are no more than 10 10 m -2 , a value that 
is at least two orders- of -magnitude less than the theoretical prediction. We 
conclude that most of the vortices detected in the original experiment must 
have been created through conventional flow processes. 

PACS numbers: 11.27. +d, 05.70.Fh, ll.10.Wx, 67.40. Vs 



1. INTRODUCTION 

After a physical system has passed rapidly through a continuous phase- 
transition, its order-parameter can have components with large differences 
between adjacent, but causally disconnected, "domains". In such systems 
topological defects can form at the domain boundaries .0 This idea was 
proposed by Kibblea in connection with the grand unified theory (GUT) 
symmetry-breakine phase-transition of the early Universe, and has been de- 
veloped by ZurekBo'Q who has estimated how the density of defects created 
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depends on the rate at which the system passes through the transition. Zurek 
also pointed out that this mechanism of defect production was applicable, 
in principle, to all continuous phase-transitions, and that it should there- 
fore be possible to validate some aspects of cosmological theories through 
laboratory-scale experiments. The first of these examined weakly first-order 
phase-transitions in liquid crystalsBQ Later, the corresponding experiments 
were carried out using thesecond-order superfluid phase-transitions of liq- 
uid 4 Hea and liquid 3 He.Bt3 All these experiments produced defect densities 
reported as being consistent with Zurek's estimates.!' BE 

In this paper we describe an improved version of the 4 He experiments 
in which particular care has been taken to minimise the production of topo- 
logical defects (i.e. quantised vortices) by ordinary hydrodynamic flow pro- 
cesses. Our new results, of which a preliminary reportiliJ has already been 
published, show no convincing evidence of any vortex creation at all. But 
they allow us to place an approximate upper-bound on the initial density of 
vortices produced by the Kibble-Zurek mechanism. 



2. THEORETICAL BACKGROUND 

A considerable amount is known about the pr op ertiesUl' 111' O^Lol vor- 
tices, and about how they are created at very low temperatures^' E-0' £2 but 
the mechanism responsible for the vorticity that appearsS'i as a result 
of passing through the A-transition (which separates the normal helium-I 
and the superfluid helium-II phases) is not understood. One possibility is 
that pre-existing rotational flow caused by e.g. convection or boiling in the 
helium-I phase is, with the onset of long-range order, converted into quan- 
tised vortex lines in helium-II; this is quite distinct from the Zurek scenario 
which we shall now briefly, describe. 

The underlying ideaB is quite simple. A small isolated volume of 
helium-I is initially held at pressure P;, and temperature TJ, just above 
the temperature T\(F\) of the A-transition. The logarithmic infinity in its 
heat capacity at T\ makes it impossible to cool the sample quickly into the 
superfluid phase but the pressure dependence of T\ means that it can be 
taken through the transition very rapidly by adiabatic expansion into the 
helium-II phase, to a final pressure Pf and temperature Tf (Fig |]). Fluctua- 
tions present in the helium-I are expected to cause the nascent superfluid to 
form with a spatially incoherent order-parameter, corresponding to a large 
density of vortex lines. This scenario depends on the fact that the liquid can, 
in principle, expand at velocities comparable to that of first-sound, whereas 
the propagation velocity for changes in the order-parameter is limited by the 
much slower velocity of second-sound. 
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The analogyBSi between liquid helium and the early Universe arises 
because they both be considered to undergo second-order phase transitions 
describable in terms of Ginsburg-Landau theory .E3 In each case the potential 
contribution to the free-energy density can be written as: 

V = a{T)\^ + (1) 

where the parameter a is positive at temperatures above T\ and negative be- 
low it, and (3 is a constant. For liquid 4 He the order-parameter is the modulus 
of the complex-scalar field ip, i.e. the Bose condensate wave-function which 
is a solution of the Ginsburg-Pitaevskii equation. In the cosmological anal- 
ogy the components of ip are Higgs fields.Bi In the symmetric (helium-I) 
phase T > T\ and the time-average of the order-parameter (tp) = 0. Be- 
low T\ this gauge symmetry, is broken so (ip) becomes non-zero, and the 
real and imaginary parts of the potential in equation 1 acquire the same 
"sombrero" shape as the corresponding cosmological free energy expressed 
in terms of Higgs fields (Fig. ^) . In the early Universe a symmetry-breaking 
phase-transition from a false- vacuum state to a true- vacuum state is thought 
to have occurred once the temperature had fallen to ~ 10 27 K, about 10 -35 s 
after the big bang. Although there are many variants of the basic modeL 
with and without inflation, it is believed that a variety of topological defectsS 
would have been produced in the transition because an event horizonnre- 
vented adjacent regions from being causally connected. Cosmic stringstil - 
thin tubes of false vacuum — are one such defect and may have had a role 
in galaxy formation. It is these that correspond to the quantised vortices 
found in helium-II. The analogy between the helium and cosmological phase 
transitions may thus be summarised as follows - 



Higgs field 1 < — > Re ip 

Higgs field 2 < ► Im tp 

False vacuum < — ► He I 

True vacuum < — > He II 

Cosmic string < > quantized vortex line 



3. DETECTING VORTICITY 

In our initial experiments the vortex density was measured by 
recording the attenuation of a sequence of second-sound pulses propagated 
through the helium-II. We expected that, following an expansion, the pulse 
amplitude would grow towards its vortex-free value as the tangle decayed 
and the attenuation decreased. 
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A considerable amount is known about the decay00lll'ii of hydrody- 
namically created vortex tangles in helium-II. Numerical simulationsBc! 
give a good qualitative description of the manner in which a homogeneous 
isotropic tangle evolves and decays. The rate at which it occurs in this 
temperature range is governed by the VinenB equation 

-TT = -X2—L 2 2 

dt TO4 

where L is the vortex-line density at time t, the 4 He atomic mass is 7714, and 
X2 is a dimensionless parameter. The relationship between vortex density 
and second-sound attenuation is knownEl from experiments with rotating 
helium, and may for present purposes be written 

L = ^ln(5„/ S ) (3) 

where C2 is the velocity of second-sound, S and So are the signal amplitudes 
with and without vortices present respectively. B is a weakly temperature 
dependent parameter, k = h/m^ is the quantum of circulation, and d is the 
transducer separation. Integrating equation ^ and substituting for L from 
equation || gives an expression for the recovery of the signal: 

I = ^2L f^t J_\ ( ) 

ln(S ) - ln(S) KBd\2iT Lj KJ 

where L\ is the vortex density immediately after the expansion. All the 
constants in this expression are known, although \2 and B do not seem to 
have been measured accurately within the temperature range of interest. 



4. THE FIRST EXPERIMENT 

A description of our first attempt to realise a bulk versionH of Zurek's 
experiment has been given in a previous paperEl but, briefly, the arrangement 
was as follows. A cell with phosphor-bronze concertina walls was filled, by 
condensing in isotopically ^1 pure 4 He through a capillary tube, and then 
sealed with a needle- valve (Fig. |||). The top of the cell was fixed rigidly to 
the cryostat but its bottom surface could be moved to compress the liquid, or 
released to expand it, using a pull-rod from the top of the cryostat. The cell 
was in vacuo surrounded by a reservoir of liquid 4 He at ~ 2K. A Straty- 
Adams capacitance gaugecll recorded the pressure in the cell and carbon- 
resistors were used as thermometers, one in the reservoir, one on the cell. 
The temperature of the cell could be adjusted by means of an electrical heater 
and a breakable thermal link to the 2 K reservoir. A trigger mechanism on 
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the mechanical linkage allowed the cell to increase its volume by ~ 20% very 
rapidly under the influence of its own internal pressure. 

As described above, the vortex density was inferred from the amplitude 
of second-sound pulses created with a thin-film heater. The signal was de- 
tected with a bolometer and passed, via a cryogenic FET preamplifier, to a 
Nicolet 1280 data processor which recorded some ~ 200 pulses during 1-2 s 
following the expansion. The last pulses in the sequence define So, the signal 
amplitude in the (virtual) absence of vortices .c3 Fig. |] shows examples of 
data recorded with this first version of the experiment. 

There was a "dead period" of about 50 ms after the mechanical shock of 
the expansion which caused vibrations that obscured the signals. The initial 
vortex-density was therefore obtained from plots such as Fig. |I](a) by back- 
extrapolation to the moment t = of traversing the transition and was found 
to be ~ 10 12 — 10 13 m~ 2 consistent with the theoretical expectation.BQH 
However, an unexpected observation in these initial experiments^ was that 
small densities of vortices were created even for expansions that occurred 
wholly in the superfluid phase (Fig. |||c), provided that the starting point 
was very close to T\. The phenomenon was initially!! attributed to vor- 
tices produced in thermal fluctuations within the critical regime, but it was 
pointed out^l that effects of this kind are only to be expected for expansions 
starting within a few microkelvin of the transition, i.e. much closer than the 
typical experimental value of a few millikelvin. The most plausible interpre- 
tation — that the vortices in question were of conventional hydrodynamic 
origin, arising from non-idealities in the design of the expansion chamber 
- was disturbing, because expansions starting above T\ traverse the same 
region. Thus some, at least, of the vortices seen in expansions through the 
transition were probably not attributable to the Zurek-Kibble mechanism 
as had been assumed. It has been of particular importance, therefore, to 
undertake a new experiment with as many as possible of the non-idealities 
in the original design eliminated or minimised. 



5. THE NEW EXPERIMENT 

An ideal experiment would avoid all fluid flow parallel to surfaces dur- 
ing the expansion. This could, in principle, be accomplished by the radial 
expansion of a spherical volume, or the axial expansion of a cylinder with 
stretchy walls. In neither of these cases would there be any relative motion 
of fluid and walls in the direction parallel to the walls, and hence no hydro- 
dynamic production of vortices. However, it is impossible to eliminate the 
effects of fluid-flow completely from any real apparatus. At the A-transition 
the critical velocity tends to zero, so any finite flow-velocity will create vor- 
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ticity. However, the period during which this happens is very short; the 
entire expansion takes only a few milliseconds so only limited growth can 
occur from the surface-vorticity sheet of half-vortex-rings.E3 

With hindsight, we can identify the_principal causes of unintended vor- 
tex creation in the original experiment]!!! in order of importance, as follows: 
(a) expansion of liquid into the cell from the filling capillary, which was closed 
by a needle valve 0.5 m up-stream; (b) expansion of liquid out of a shorter 
capillary connecting the cell to the Straty- Adams capacitive pressure-gauge; 
(c) flow out of and past the fixed yoke (a U-shaped structure) on which the 
second-sound transducers were mounted. In addition, (d) there were the 
undesirable transients caused by the expansion system bouncing against the 
mechanical stop, at room temperature. The walls of the cell were made from 
bronze bellows,&B3 rather than being a stretchy cylinder, but the effects of 
the flow parallel to the convoluted surfaces were relatively small. 

The design of our new expansion cell (Fig. addressed all the listed 
problems, as follows: (a) a hydraulically-operated needle-valve eliminated 
the dead- volume of capillary tube; (b) the phosphor-bronze diaphragm of the 
pressure gauge!! became an integral part of the upper cell-wall eliminating 
the long tube leading to the pressure gauge; (c) the cell was shortened from 
25mm to ~ 5mm so that the heater-bolometer pair could be mounted flush 
with its top and bottom inner surfaces, thereby eliminating any need for a 
yoke; (d) some damping of the expansion was provided by the addition of a 
(light motor-vehicle) hydraulic shock-absorber. 

The operation of the apparatus, the techmcme of data collection and 
the analysis were much as described previously^' except that the rate at 
which the sample passed through the A-transition was determined directly 
by simultaneous measurements of the positionlll of the pull-rod (giving the 
volume of the cell, and hence its pressure) and the temperature of the cell. 

The position detector made use of the variation of inductance of a 
solenoid coil when a ferrite-core is moved in and out. The 1 cm long, 4 
mm diameter, ferrite-core was rigidly attached to the top of the pull-rod, 
where it emerged from the cryostat, and was positioned so that half of it 
was inside the 5 mm bore of the solenoid. The solenoid itself was clamped to 
the cryostat top-plate so that, as the cell expanded and the pull-rod moved 
down, the ferrite penetrated further inside. The position was measured by 
making the solenoid part of a resonant LCR circuit and setting the frequency 
so that the output was at half maximum. The 2 MHz AC signal was rectified, 
smoothed and then measured using a digital oscilloscope. The response time, 
which was limited by the smoothing and the Q of the circuit, was ~ 50/US 
and the sensitivity was ~500 mV/mm. This, and the linearity of the system, 
were measured by calibration against a micrometer. The position could in 
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principle be measured to an accuracy of ~ ±2/j,m but in the experiment this 
became ~ ±Wfim because of the 8 bit digitiser of the oscilloscope. 
We define the dimensionless distance from the transition by 

T x -T 



(5) 



and the quench time r q by 



-) 



) (6) 

e=0 



Fig. ^ shows a typical evolution of e with time during an expansion. In 
Fig. ^(a), which plots the full expansion period, it is evident that the system 
"bounces" momentarily near e = 0.02, but without passing back through the 
transition. The effect is believed to be associated with the onset of damping 
from the shock absorber, after backlash in the system has been taken up. 

The quench time is readily determined from e(t) near the transition. In 
the case illustrated in Fig. [f| shown in expanded form in part (b), it was 
r q = 17 ± 1 ms. We are thus able take both the pressure-dependence of T\, 
and the non-constant rate of expansion, explicitly into account. 



6. RESULTS FROM THE NEW CELL 

The improvements in the cell-design intended to eliminate vortices orig- 
inating from conventional fluid-flow were clearly successful; as we had hoped, 
and unlike the first cell, expansion trajectories that stay within the super- 
fluid phase, even those starting as close as 42 mK to the transition, create 
no detectable vorticity (Fig. 0). 

We were surprised, however, to discover that no detectable vorticity was 
created even when the expansion trajectory passed through the transition. 
Signal amplitudes measured just after two such expansions are shown by 
the data points of Fig. || It is immediately evident that, unlike the results 
obtained from the original cellja there is now no evidence of any systematic 
growth of the signal amplitude with time or, correspondingly, for the creation 
of any vortices at the transition. One possible reason is that the density of 
vortices created is smaller than the theoretical estimates,!' i'i but we must 
also consider the possibility that they are decaying faster than they can be 
measured. 



7. THE DECAY OF A VORTEX TANGLE 



To try to clarify matters, we performed a subsidiary set of experiments, 
deliberately creating vortices by conventional means and then following their 
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decay by measurements of the recovery of the second-sound signal amplitude. 
By leaving the needle-valve open, so that ~ 0.2 cm 3 of liquid from the dead 
volume outside the needle- valve actuator-bellows squirted into the cell during 
an expansion, we could create large densities of vorticity and observe their 
decay. Despite the highly non-equilibrium situation that arises immediately 
following the expansion, as liquid squirts into the cell, it was found that the 
temperature reached a steady value after ~ 6 ms. Fig. |9] shows two examples 
of attenuation plots ([ln^/So)] -1 against time t) and these always, within 
the experimental errors, had the linear form predicted by equation [|. Fig. [K| 
summarises the results of a number of expansions from which it was possible 
to determine X2/B as a function of temperature and pressure. We found 
that X2/B was weakly temperature-dependent and over the range of interest, 
0.02 < e < 0.06 it could be approximated by X2/B = 0.004 ± 0.001. 
This measured value of X2/B was then used to calculate the evolution of 
S/Sq with time for different values of Li, assuming B = 1, yielding the 
curves shown in Fig. 8. From the r q derived from the gradient in Fig |6|, and 
Zurek's estimate (based on renormalisation-group theory) of 

L RG = L ° where L = 1.2 x 10 12 m~ 2 , r = 100 ms (7) 

we are thus led to expect that L\ ~ 4 x 10 12 m 2 . Direct comparison of the 



calculated curves and measured data in Fig. 11 shows that this is plainly 
not the case. In fact, the data suggest that L\, the vorticity created by the 
transition, is no more than 10 10 m~ 2 , smaller than the expected value by at 
least two orders-of-magnitude. 



8. DISCUSSION 

Given the apparently positive outcome of the earlier investigations,! 
the null result of the present experiment has come as something of a sur- 
prise. There are several points to be made. First, Zurek did not expect his 
estimates of L\ be accurate to better than one, or possibly two, orders-of- 
magnitude, and his more recent estimatell suggestesiirnewhat lower defect- 
densities. So it remains possible that his pictured' qBl3 is correct for 4 He in 
all essential details, and that an improved experiment with faster expansions 
now being planned will reveal evidence of the Kibble-Zurek mechanism at 
work in liquid 4 He. Secondly, it must be borne in mind that (g), and the 
value of x/B measured (from plots like those in Fig. [n]) from the data of 
Fig. ^, refer to hydrodynamically generated vortex lines. Vorticity generated 
in the nonequilibrium phase transition might perhaps be significantly differ- 
ent, e.g. in respect of its loop-size distribution.! 33 ! It could therefore decay 
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faster, and might consequently be unobservable in the present experiments. 
Thirdly, it is surprising to us that the 3 He experiments seem to agree 
with Zurek's original estimatesBcl'El whereas the present experiment shows 
that they overstate L\ by at least two orders-of-magnitude. It is not yet 
known for sure why this should be, although an interesting explanation of 
the apparent discrepancy has recently been suggested^ by Karra and Rivers. 
They suggest that fluctuations near to T\ may change the winding number, 
i.e. reduce the density of vortices produced in the 4 He experiment, and they 
show that the analogous density reduction would be much smaller in the 
case of the 3 He experiments. 
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Fig. 1. Schematic of expansion trajectory through the 4 He superfluid transi- 
tion from a starting temperature and pressure (Tj, Pj) to final values (T* , PA. 
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Fig. 3. Schematic diagram showing the main features of the original exper- 
imental cell. The heater and bolometer were mounted on a yoke immersed 
in the liquid 4 He. 
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Fig. 4. Recovery of the scaled second-sound signal amplitude S as a function 
of time t in our original experiment, for various expansion trajectories: (a) 
through the lambda line, 7] = 1.81 K, P; = 29.6 bar, T f = 2.05 K, P f = 
6.9 bar; (b) starting well below the lambda line, T\ = 1.58 K, P = 23.0 bar, 
T f = 1.74 K, P f = 4.0 bar; (c) starting slightly (~ 10 mK) below the lambda 
line, 7] = 1.82 K, P; = 25.7 bar, T t = 2.03 K, P f = 6.9 bar. 
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valve 



Fig. 5. The new expansion cell, designed so as to minimise hydrodynamic 
creation of vortices. 
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Fig. 6. Distance versus time t for a typical quench: (a) during the complete 
expansion; (b) enlargement of region near e = 0. The reciprocal of the 
gradient at the transition gives the value of the 'quench time' parameter, r q . 
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Fig. 7. Second-sound pulse amplitude S/ Sq as a function of time following 
pressure quenches entirely within the helium-II phase: (a) starting far (~ 
490 mK) below T x with T { = 1.37 K, P ; = 24.2 bar, T f = 1.47 K, P f = 7.2 bar; 
(b) starting ~ 150 mK below T A with T { = 1.74 K, P { = 23.8 bar, T t = 1.76 K, 
P f = 6.9 bar; (c) starting slightly (~ 40 mK) below T\ with Tj = 1.85 K, 
Pi = 22.4 bar, T t = 1.98 K, P t = 6.4 bar. 
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Fig. 8. Second-sound pulse amplitude S/So as a function of time follow- 
ing pressure quenches through the A transition: (a) starting close to the 
transition, with T; = 1.81 K, P = 30.3 bar, T f = 2.03 K, P f = 6.2 bar 
(ej = —0.032); (b) starting further above the transition, with 7] = 2.05 K, 
Pi = 22.7 bar, T f = 2.09 K, P f = 5.9 bar (e; = -0.089); (c) starting far above 
the transition, with T| = 1.96 K, P; = 34.1 bar, T f = 2.07 K, P f = 6.1 bar 
( €i = -0.167). 
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Fig. 9. Examples of (lnS - lnS")" 1 plotted versus time t for hydrodynami- 
cally created vortices. 
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Fig. 10. Measured values of the parameter X2/-B as a function of distance e 
from the A-transition. 
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Fig. 11. Evolution of the second-sound amplitude S with time, following an 
expansion of the cell at t = (data points), normalised by its vortex- free 
value Sq. The open and closed symbols correspond to signal repetition rates 
of 100 and 50 Hz respectively, and in each case the starting and finishing 
conditions were = — 0.032, e/ = 0.039. The curves refer to calculated signal 
evolutions for different initial line densities, from the bottom, of 10 12 ,10 11 
and 10 10 m~ 2 . 



